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Abstract 

The general form of the stress-tensor three-point function in four dimensions is 
obtained by solving the Ward identities for the diffeomorphism and Weyl symme- 
tries. Several properties of this correlator are discussed, such as the renormalization 
and scheme independence and the analogies with the anomalous chiral triangle. At 
the critical point, the coefficients a and c of the four-dimcnsional trace anomaly are 
related to two finite, scheme-independent amplitudes of the three-point function. Off- 
criticality, the imaginary parts of these amplitudes satisfy sum rules which express the 
total renormalization-group flow of a and c between pairs of critical points. Although 
these sum rules are similar to that satisfied by the two-dimensional central charge, the 
monotonicity of the flow, i.e. the four-dimensional analogue of the c-theorem, remains 
to be proven. 
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1 Introduction 



The Zamolodchikov c-theorem [|| states that the two-dimensional unitary field theories 
are characterized by a positive function of the coupling constants c{g) which decreases 
along the renormalization-group flow; moreover, its value at the fixed points is equal to 
the trace-anomaly number c. Several works have suggested that the theorem could extend 
to four dimensions^: at criticality, the candidate c-function would again reduce to the 
coefficient of the the Euler density in the trace anomaly |p. This conjecture is supported 
by many examples of non-trivial flows in supersymmetric gauge theories, which have been 
analysed via the Seiberg-Witten electric-magnetic duality and the AdS supergravity- 
gauge theory correspondence |8| . Two proofs of the four-dimensional theorem have also 
been proposed, although one is rather controversial and the other one relies on a strong 



dynamical hypothesis |10|. 



These results and conjectures motivate further analyses of the physical consequences 
of the four-dimensional trace anomaly^ 

g^'^iT^u) ^{Q) = X{aE-3cW + a' D'n + rTZ') , A ^ ^^^^^ • (1-1) 

In this equation, E denotes the Euler density (x oc / E is the Euler characteristic), 
W the square of the Weyl tensor and TZ the curvature scalar (see Appendix ^ for our 
conventions). Some properties of the trace anomaly ( |1.1| ) are known since a long time 
[|^][|^]: the fact that the Wess-Zumino consistency conditions imply r = |14]|15| and 



that the coefficient a' is scheme dependent, because it can be changed by adding a local 



counterterm to the effective action [13|. An interesting study of the four-dimensional trace 



anomaly (and its higher-dimensional analogues) has been done in Ref.|16|, by deriving 
the corresponding terms in the regularized effective action within the framework of the 
epsilon-expansion. It was found that the topological invariant Euler term corresponds to 
an effective action which is finite and non-local as e ^ 0; this is the same feature of the 
chiral anomaly |12|. On the other hand, the W term in (1.1) corresponds to a singular 



effective action with finite scale variation. 

In this work, we relate the four-dimensional trace anomaly to some physical quantities: 
i) we express the universal parameters a and c in terms of finite scheme- independent 
correlation functions at criticality; ii) away from criticality, we provide exact sum rules for 
the renormalization-group flow[| of a and c. 

* For a short review, see Ref.Q; for early attempts, see the Refs.|^] 

^ Note that the normahzation of the four-dhnensional stress tensor in Eq.(l.l) differs from that in 
Refs.[^[0 by the factor l/47r^; nevertheless, the values of the anomaly coefficients for the free conformal- 
invariant scalar field are still normalized to a = c = 1. 

* These flows have also been analysed by the Refs. |4| |tLO| within the perturbative epsilon-expansion. 
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The analogous two-dimensional results are well known ll^]: the coefficient c of the 
trace anomaly (0) = — c 7^/12 is unambiguously defined by the stress-tensor two-point 
function (T(z)r(O)) = c/2z^ at criticality, which is a finite and well-defined correlator. 
Furthermore, the renormalization-group flow of c is constrained by the Cardy sum rule, 



which involves the off-critical correlator |18]: 



3 I' f°° 
CUV -ciR = — / d^x x'^{Q{x)Q{0)) = / ds c{s) > . (1.2) 

J\x\>e Je 

In the above Equation, cjjv and cjr are the central charges of, respectively, the ultraviolet 
and infrared fixed points connected by the renormalization-group trajectory, and the x- 
integral is evaluated at any point along the trajectory. Equation ( |1.2D also shows the form 
of the sum rule in momentum space H]: it involves the positive-definite spectral measure 
c{s) ds, which is obtained from the imaginary part of {Q{p)@{—p)) at momentum = s. 
In four dimensions, the relation of the trace anomaly ( |1 . 1| ) to finite correlators requires 
the study of the three-point function in flat space|^ or of the two-point function in curved 
space; their general expressions at criticality (in coordinate space) have been found in the 
Refs.|l|l@ and 0, respectively. 



In this paper, we obtain the general form of the three-point function in (flat) momentum 
space at and off criticality: we write it as a sum of scalar amplitudes (form factors) 
times tensor structures, which solve the Ward identities for the diffeomorphism and Weyl 
symmetries (the latter being imposed at criticality). In addition, we write dispersion 
relations for the renormalized amplitudes following closely the well-known analysis of the 
chiral triangle [p2|. 



At criticality, we can identify the three-point amplitudes which correspond to the terms 
in the trace anomaly, disentangling the scheme-dependent amplitude related to a' in ( p..l| ) 
from the scheme-independent ones corresponding to the coefficients a and c. Any ampli- 
tude can be singled out from the general expansion of (TTT) by solving a linear system: 
we obtain some reduction formulae which can be useful for computing a and c in interact- 
ing theories; these relations could also suggest further physical consequences of the trace 
anomaly. 

Away from criticality, the imaginary parts of the a and c amplitudes satisfy two exact 
sum rules, which express the total variation of a and c under the renormalization-group 
flow between two fixed points; these sum rules follow from the analysis of the general 
structure of (TTT) and are analogous of the two-dimensional result ( |1.2[ ). 

In two dimensions, the imaginary part of the two-point function c(s) is positive definite 



by unitarity; the sum rule (1.2) then implies that c decreases along the renormalization- 



Hereafter, we shall mostly use the Euclidean signature (see Appendix A for our conventions). 
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group flow. In four dimensions, tlie property ajjv > ajR could be similarly proved if 
the imaginary part of the a amplitude were positive definite. In general, we do not 
expect that three-point functions, related to decay amplitudes, should have a definite sign; 



nevertheless, we should mention the positivity condition proposed in Ref.[20|, coming from 
the classical weak-energy condition [p^j, which implies the positivity of the a coefficient at 
criticality. The consequences of this condition on the three-point amplitudes off-criticality 
remain to be analysed. 

The plan of the paper is the following: in Section 2, we briefly recall the dispersive anal- 
ysis of the chiral triangle (AW) of Ref . [p2|| ; this lets us outline some features which arise 
for (TTT) too. In this Section, we also give the definition of the stress-tensor operator and 
repeat the dispersive analysis of (TT) ||5|[24|p5|. Appendix ^ and |^ contain our notations 



for the curved-space calculus and the derivation of the Ward identities, respectively. 

In Section 3, we present the analysis of (TTT): i) the decomposition in scalar am- 
plitudes times the tensors solving the diffeomorphism Ward identity (the tensor basis is 
constructed in Appendix |^) ; ii) the conditions on the amplitudes which are imposed by 
the (anomalous) Weyl symmetry at criticality. Owing to the large number of possible 
tensor structures, we make an extensive use of algebraic programs (longer formulas are 
given in Appendix ^ and the full output is available as a Mathematica |26] notebook). 
In this Section, we also discuss the Ward identities that follow from flat-space conformal 
invariance at criticality. 

In Section 4, we show how to extract the scalar amplitudes in the (TTT) expansion from 
an explicit expressions of this correlator; in particular, we obtain the reduction formulae 
which determine the trace anomaly coefficients a and c at criticality. Next, we discuss 
the renormalization-group flow of the corresponding amplitudes off-criticality and derive 
the sum rules analogues to Eq.( |1.2| ). In Section 5, we test the reduction formulas and 
the sum rules by computing the flows of the free massive scalar and fermion theories. In 
the Conclusions, we briefly discuss the possible attempts towards the c-theorem in four 
dimensions. 

2 Dispersive analysis of the chiral triangle and of the stress- 
tensor two-point function 
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2.1 Analysis of {AW) 



We start by recalling the well-known analysis of the three-point function of one axial and 
two vector currents p^ : 

We shall discuss the renormalized correlator, i.e. we shall not explicitly introduce bare 
expressions and counterterms, although some finite scheme-dependent terms will occur. 
The space-time symmetries allow to decompose ( p.lj ) in tensor structures times scalar 
amplitudes Ai = Ai{q'^, fcf , fcl); we consider the symmetric point kf = in order to fully 
exploit the Bose symmetry (k"^ <^ ^2^) of the correlator. The tensor notation can be 
simplified by contracting the Lorentz indices with general polarization vectors {v^ , e^, Cg}, 



following Ref.[27|. We then write: 



r(t), 61,62) = e"' 62^ r^aiaa = ^1 bei62(A;i - k2)] 

+ A3 {[veikik2] {e2.ki) + [ve2k2ki] {ei.k2)) 

+ A5 {[veikik2] (62.^2) + [ve2A;2A:i] (ei.fci)) 

+ A7[eie2kik2]{v.{ki + k2)) , (2.2) 



ttf^b'' = (a.b) , e^yp^a^'^cPd" = [abed] . (2.3) 



where we denoted: 



The tensor structures in ( |2.2| ) are not linearly independent in four dimensions, due to the 
identity 



(a./) [bcde] + (6./) [cdea] + (c./) [deab] + {d.f) [eabc] + (e./) [abcd\ = , (2.4) 
which holds for any six-plet of vectors {a, 6, c, d, 6, /}. In its Bose symmetrized form, this 



identity implies that one of the scalar amplitudes in the expansion (2.2) is redundant. 



Actually, T is invariant under the shift of amplitudes by an arbitrary function /: 

Ai^Ai- {e + {kiM)) / , Ai^Ai-f, for i = 3, 5, 7 . (2.5) 



Therefore, we can eliminate one of the amplitudes in (2.2): we choose = to match 



the notations of Ref.|22 



We consider a general theory with massive fermions and write the Ward identities for 
the conservation of the vector and axial currents which enter in {AW) (the axial current is 
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supposed anomalous). These Ward identities can be expressed in our notation by choosing 
longitudinal polarizations in T (v, 61,62), for example ei = ki. We then write, 



T{v,ki,e2) = , 

(J 

r (fci + /c2, ei, 62) = 2m Ts (61,62) + — [6162/^1^2] ; (2.6) 

vr 

^5 (^1,62) is the Fourier transform of the corresponding triangle (J^VV), involving the 
pseudo-scalar density, whose tensor structure is: 

Ts (61,62) = - vlg [6162^1/02] . (2.7) 

The resulting conditions on the scalar amplitudes are: 

2yli = g2^3 + 2A;2 (^5 - A3) ; (2.8) 
C 

2Ai = h2m Aq . (2.9) 

TT 



The quantity C on the right-hand side of ( |2.9D is the chiral anomaly: C cannot be put 
to zero because T would vanish altogether for m = at the point q'^ = k'^ = 0, which 
is not physically acceptable. Furthermore, C can be shown to be a constant as follows 
[p^ ]. Consider the imaginary parts of all amplitudes w.r.t. both variables and k"^: they 
are finite quantities describing physical scattering processes, that satisfy non-anomalous 
Ward identities, 

2ImAi(g^A;2) = (g^ - 2A;2) Im ^3(^2, fe^) + 2A;2 Im ^(g^ A;^) , (2.10) 
ImAi(g2,A;2) = mlm A9{q^,k'^) . (2.11) 

It follows that the anomalous term has vanishing imaginary parts: 

C(g2 + ie, k^) - C{q^ - ie, k'^) = = C{q^, k^ + ie) - C{q^, k^ - ie); (2.12) 

thus, C is a dimensionless polynomial of /m? and kP' jrr? with finite m — > limit, i.e. it 
is a constant. 

Next, we write the dispersion relations w.r.t. g^: 

AM,k-)=A\^^^-\,s ^-^Msf) ^ (2.13) 
■K J s-q^ 

AM\k^) = - jds ^'^^^^'f) , i = 3,5,9. (2.14) 
TT J s-q^ 

The integrals in Eq. (2. IS ,2.14) are convergent due to the asymptotic behaviors, 

Ai{q^, k^) ~ q '^''^A' , for q^ ^ 00 , (2.15) 
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following from dimensional analysis (up to logarithms which are supposed absent for ^i): 
dimAi = 0, dimAa = dimAs = —2 and dimAg = —1; only the amplitude Ai (superficially 
divergent at the bare level) requires the subtraction term A^^ in ( 2. 131 ). Combining Eqs. 
(|1|,|1|) and (|1|), we can derive the sum rule: 



^Jdslm As{s, fc2) + (g2 _ 2A:2) ^3(^2, + 2A;2^5(g2, fc2) = 2mAg{q^, /fc^) . (2.16) 



The comparison of this sum rule with the original Ward identities ( p.^j2.9| ) determines 
Ai^ in ( 2.13| ) as the chiral anomaly C and yields the following sum rule: 

= = Ms,k') , (2.17) 



Let us further analyse the Eqs. (2.8, 2.9 , 2.17| ) for k'^ = m? = 0: A^ and A^ drop out 
and we find the explicit solution, 

A^{q\^) = lmA^{q\Q)=C6{q^), 

^i(g',0) = -^ = -l-J dslmA3is,0). (2.18) 



This is famous result ||2^ that the anomaly appears as a delta term in the imaginary 

part of a finite amplitude (^3). We also notice that in the general case 

the constancy of the sum rule ( 2.171) implies a non-trivial constraint on Im A^{q'^, k'^,m): 



this should be a smoothed delta-function with spreading parameter proportional to k"^ for 
m? = (respectively, m? for k"^ = Q). 

Let us summarize a number of facts which will also occur in the analysis of (TTT): 

• There could be algebraic degeneracies among the tensor structures which can be 
written for the three-point function. 

• The anomaly can be expressed in terms of a finite, scheme-independent amplitude 
of the three-point function (i.e. A3{q'^,0) oc C/q^); this relation is analogous to 
{T(z)T(0)) = cl2z^ in two dimensions. 

• Scheme-dependent amplitudes (i.e. A\) are related to scheme-independent, finite 
amplitudes (i.e. ^3) by the Ward identities. 



The sum rule for the chiral anomaly (2.17) is similar to the sum rule (|1.2| ) for 



the conformal anomaly in two dimensions, which is a version of the c-theorem: both 
expressions involve dimensionless densities which are smoothed delta- functions. This 
similarity will be fully explained in Section 4, where the sum rules for the four- 
dimensional trace anomalies will be derived by combining both arguments. 
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2.2 Definition of the stress tensor and Ward identities 



Here we consider a four-dimensional renormalizable field theory and suppose that it can 
be covariantly extended on classical curved-space backgrounds; at critical points, we also 
assume that the theory is covariant under Weyl transformations of the background metric 
9fiu{x)- The stress-tensor n-point function is defined by the following variation of the 
generating functional W[g] = log Z[g] w.r.t. the metric: 

1 6 6 

{T,,,', (xi) • • • r,„,,„ ^ ^^...^/^ 6g>^^^'.{x^) ■ ■ ■ 6g^->^'^{x„)^^^^ ' 

(2.19) 

In this Equation, the ^ terms are pulled out of the variations, so that n-point functions 
explicitly respect the symmetry under exchanges of two stress tensors. Other definitions 
would yield n-point functions differing by lower /c-point functions, < < n; this is an 
intrinsic ambiguity of the stress-tensor which should be fixed together with the renormal- 
ization conditions (this point will be further discussed in Section 5). 

The symmetries of coordinate reparametrizations (Diff) and, at critical points, of Weyl 
rescalings of the metric imply Ward identities which can be obtained by varying W[g] 
and by consistently using the definition of the stress tensor ( p. 19 ) (see Appendix ^ for 



the conventions of curved-space calculus and Appendix ^ for a complete discussion of the 
Ward identities). We can write: 

Diff : 6,W[g] ^ [ dx {D^'e^x) + D^e^ix)) —^W[g] = ; 

J 6g^"'{x) 

Weyl : 6„W[g] ^ - I dx 2a{x) g^^- {x)—^W[g] 

J 6gt"'{x) 

= - jdx 2aix) ^Aix;g) . (2.20) 

Note that the Weyl symmetry is anomalous with trace anomaly A{x; g). The first variation 
w.r.t. g^^y of the identities ( p. 201) yield the Equations: 

Diff: D^'{T^,,>)=Q; (2.21) 
Weyl: g^^' {T^^,) = {Q) = A{x;g) = \{a E - 2>c W + a' D'^U + r U^) . (2.22) 

The Ward identities for multi-point functions are obtained by further variations w.r.t. 
g^^'{x) of the Eqs. (g^Mj) ([l22|) . 



The expression ( ^.221) of the trace does not contain explicit dynamical fields which 
would violate the Weyl symmetry at the classical level {Q classical = 0): the underlying 
hypothesis is that the scale-invariant theory in flat space can be extended to a Weyl- 
invariant theory in curved space. Above two dimensions, this extension is not automatic: 
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flat-space scale invariance is compatible with derivative fields in the trace, which, however, 
can often be removed by redefining (i.e. "improving" ) the stress tensor [^] (obstructions 
to improvement are discussed, e.g., in Ref.p). 



2.3 The form of (TT) and its scheme dependencies 

The general form of the stress-tensor two-point function is one of the ingredients which 
are needed for the construction of the three-point function: we recall here the analysis of 
Ref. Q, which can also be found in |13]|25]. From the variation of ( 2.21| ), we obtain the 



diffeomorphism Ward identity for the two-point function in flat space: 
2^^{T,,iy)T^p{x)) + ^ {5iy - x)(r^),(x))) + (^<5(y - x)) (r„^(x)) = , (2.23) 

where the symmetrization of indices is x^^y'^^ = x^y'^ + y^x'^ . The standard normalization 
condition would fix the constant {T^u) = in flat space; at critical points, this is consistent 
with the absence of dimensionful parameters in the trace anomaly, such as a cosmological 
constant. On the other hand, we prefer to keep a non- vanishing one-point function for 
the later purpose of studying the scheme-dependence of (TTT) due to off-critical mass 
parameters; thus, (T^jjy) oc rji^ty m^. 

The first variation of ( ^.22 ) similarly yields the Ward identity for the Weyl symmetry 



at criticality: 

(G(y)T«^(x)) = Xa' (r/a/jS^ - do^dp) d'^5{x - y) , (critical points). (2.24) 

It is convenient to contract the indices of the stress-tensor with generic polarization tensors 
h^u = hy^, and use the following notations: 

K = (^) ' {hi)ap{h2r^ = {hiM) , v'^Kpq^ ^ (p.h.q) . (2.25) 

The conservation and trace equations in momentum space will be obtained by respectively 
choosing the polarizations: 

hi^^{v<S)ki + ki<S)v) , hi^r], « = 1,2, (2.26) 

where v is an arbitrary four- vector, ki is the momentum of the Fourier transformed tensor 
T{ki) and ry the Euclidean (Minkowskian) metric (our conventions are listed in Appendix 

0)- 

In conclusion, the Diff and Weyl Ward identities for the two-point function in momen- 
tum space read: 

Diff : {{v p.T{p)) {h2.T{-p))) + 2{{p.h2.T.v)) + {v.p){{h2.T)) = , (2.27) 
Weyl: {Q{p) {h2.T2{-p))) = \ a' p^ {p^{h2) - {pM-p)) , (critical points). (2.28) 
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The general form of the two-point function in momentum space is obtained by the 
procedure of Section 2.1: the general expansion in tensor structures times scalar amplitudes 
is constrained by the Ward identities. The one-point function can be written: 

{(h.T)) = -2 {h) /a , (2.29) 

where f\ is a dimensionful constant related to the mass. The two-point function solution 



of the Diff Ward identity (|2.27D is found to be l5|: 



{{h,.T{p)) {h2.n-p))) 

= (/O {p') - f2 {p\h) - {p.h,.p)) {p\h2) - {p.h2.p)) 

+ (d- l)/2 (p^) [{pM-P){V-h2.p) - 2p\p.hi.h2.p) +/(/il./l2)] 

+ /a [(/ii)(/i2) + 2(/ii./i2)] . (2.30) 

The two independent scalar amplitudes /o (p^) and /2 (p^) are dimensionless functions of 
the squared momentum and the mass; they parameterize the propagation of intermediate 
states of spin zero and two, respectively, and are multiplied by transverse tensors. More- 
over, the spin- two tensor is traceless (in two dimensions, it vanishes identically providing 
another example of tensor degeneracy). For later purposes, it is convenient to write the 
solution for general space-time dimension d, and set d = 4 when comparing with the trace 
anomaly. 

At the critical points (m = 0), the two dimensionless amplitudes can depend logarithmi- 
cally on the momentum or be constant: the spin-zero amplitude is completely determined 



by the Weyl Ward identity (p^ : 



[p^) = , (critical points). (2.31) 

The absence of a logarithmic term is due to our starting hypothesis of Weyl invariance 
at criticality, which forces the trace for the stress tensor to vanish as an operator; its 
correlators can only contain scheme-dependent contact terms. Thus, the spin-zero part 
of (TT) is a polynomial of momenta of dimension four which can be modified by adding 
counterterms to the effective action. This verifies the well-known scheme dependence 
a' — > a' -|- const in the trace anomaly (LI) Wi 



The amplitude /2 (p^) has instead a logarithmic singularity at the bare level; its renor- 
malized expression is: 

/I (^2) = -A c log , (critical points). (2.32) 

In the above Equation, /x is the renormalization scale and c is the coefficient of the Weyl 
term in the trace anomaly ( p. 22 ): the latter fact can be proved ||l^||5| by using the Ward 
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identity for scale invariance in flat space (a conformal isometry fully discussed in Section 
3.3). Equation ( p. 32 ) shows that /2 (p^) is defined up to a scheme-dependent constant, 
while the coefficient c is scheme- independent: its value normalizes the two-point function 
at criticality (for non-coincident coordinates, xi ^ X2) and is positive for unitary theories. 

In conclusion, the two (TT) amplitudes depend on the renormalization scheme by 
constant shifts, which change the correlator by polynomials of the momenta. The scheme- 
dependence of any correlator can be characterized as follows: renormalization theory says 
that counterterms in the action are polynomial of momenta with maximal degree d, thus 
scheme-dependent terms of a A-dimensional n-point function {Oi{ki) ■ ■ ■ On{kn)) are poly- 
nomials of maximal dimension A (A > is assumed here). These polynomials can be 
added to the n-point function, which is being expanded in tensor structures times scalar 
amplitudes, as done for {AVV) and {TT); if a polynomial matches one of the tensors in 
the expansion, the corresponding scalar amplitude is redefined. 

This scheme dependence can be characterized by a number associated to the corre- 
sponding tensor structure, which will be called the tensorial dimension 6t- This is de- 
fined as the highest number of free-index momenta of the monomials contained in the 
tensor (as usual, free indices are contracted with polarizations hi or v). For instance, 
5T[{ki.h2.k2)ik2.hi.k2) + k'^q'^{hi.h2)] = 4 and 6T[k'^q'^{hi.h2)] = 0. Therefore, the /o,/2 
amplitudes in (|2.3[1| ) are both associated to 5^ = 4 tensors, while /a multiplies a 6t = 
expression. 

A sufficient condition for an amplitude to be scheme independent is that its tensor 
has dimension 6t > A: clearly, tensors with (5t > A cannot mix with counterterms, 
while (5t < A tensors could mix with them - a detailed analysis is then necessary. For 
the n-point functions of the stress tensor (A = d), the scheme- independence condition is 
5r > 4 in four dimensions; all amplitudes in (TT) have associated tensors with 5^ < 4 
and are actually scheme-dependent. On the contrary, (TTT) can have 6t = Q tensors and 
corresponding finite, scheme-independent amplitudes; finding the finite amplitudes related 
to the anomaly coefficients a and c is the main motivation for studying the three-point 
function. 



3 The stress-tensor three-point function 
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3.1 Derivation of the general form 



In this Section, we derive the general form of the stress-tensor three-point function at the 
symmetric point k"^ = k\ = k"^: 

{{h^.Tik^)) {h2.T{k2)) {hi.T{ki))) , q^' = -k^ = k>l + k^ , kl = kl = e . (3.1) 

Although the point q"^ = k'^ is the most symmetric, it will be not consider here because 
it lays in the unphysical region of complex momenta, while the choice q^ 7^ k'^ allows for 
real (Minkowskian) momenta both at and off criticality. 

The diffeomorphism Ward identity constraints the longitudinal component of the stress 
tensors: for T^k^), it reads (see Appendix]^ for the derivation): 

= 2{{k3.T{k^).v){T{k2).h2){T{ki).hi)) 

+ 2{{k^.h2.T{-ki).v){T{ki)M)) - {k2.v){{T{-ki).h2){T{ki).hi)) 

+ 2{{k^M.T{-k2).v){T{k2).h2)) - {ki.v){{T{-k2)M){T{k2).h2)) . (3.2) 

Note that this inhomogeneous equation relates the three-point function to the two-point 
correlator already found in Section 2.3. 

The general form of (TTT) off criticality can be expanded over the basis of the Lorentz 
tensors which are built out of kf , k2 , rj'^'^ : there are 137 elementary polynomials Vi = 
Vi{kf, k2,'r]^'^) of maximal tensorial dimension 6t = 6, which are classified in Appendix^. 
After imposing the Bose symmetry (hi, ki) (/12, k2) at the symmetric point kf = k2, the 
number of independent polynomials reduces to 77 (see Table [l|). This general expansion is 
plugged into the Ward identity (3^) and the corresponding one for (ki ^ k^), which yield 
linear systems for the scalar amplitudes (the algebra is done with the help of Mathematica 
[ p^ ] routines). In order to avoid the problem of possible degeneracies of the tensor basis due 
to relations proper to four-dimensions, we solve the system in arbitrary large dimensions, 
and postpone the discussion of the limit to four dimension. The solution of the Ward 
identity (3^) and its (1 ^ 3) analogue contains 17 transverse tensors, solutions of the 
homogeneous equations, plus 5 inhomogeneous (i.e. non transverse) tensors multiplying 
the two- and one-point amplitudes {/o (g^) , /o (fc^) , /2 {q^) , /2 {k'^) , /a} already found in 
Section 2.3: 

17 

{{hi.T){h2.T){h^.T)) = Y.A,{q^k')% 

i=l 

+ fo{q') Tog + fo{k^) %k + f2{q^) T2, + f2 {k') 

+ /a Ta . (3.3) 

Note that the tensor structures are contracted with the polarizations (/ii, /i2, /is), i.e. 
= %{ki,hi,k2, /i2, hs), for w = 1, . . . , 17,0q,0k, 2q, 2k, A. 
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N 


N+1 


N+2 


N+3 


N+4 





(1|1|1)(1|2|1)(1|3|1) 


(1|3|1)(2|1|2)(2|2|2) 


(1|1|1)(1|2|1)(1|3|2) 


(1|1|1)(1|3|1)(2|2|2) 


4 


(1|1|1)(1|3|2)(2|2|2) 


(1|1|1)(1|3|1)(1|2|2) 


(1|3|1)(1|1|2)(2|2|2) 


(1|1|1)(1|2|2)(1|3|2) 


8 


(1|2|1)(1|3|1)(2|1|2) 


(1|2|1)(1|3|2)(2|1|2) 


(1|3|1)(1|2|2)(2|1|2) 


(1|2|1)(1|3|1)(1|1|2) 


12 


(1|2|1)(1|1|2)(1|3|2) 


(1|3|1)(1|1|2)(1|2|2) 


(1|1|2)(1|2|2)(1|3|2) 


(1|1|1)(1|32|1) 


16 


(1|31|1)(2|2|2) 


(1|1|1) (1|32|2) 


(1|1|1)(2|32|1) 


(1|32|1)(2|1|2) 


20 


(1|2|1)(1|31|1) 


(1|2|1) (2|31|1) 


(1|2|1)(1|31|2) 


(1|32|1)(1|1|2) 


24 


(1|31|1)(1|2|2) 


(1|1|2) (1|32|2) 


(1|2|2)(1|31|2) 


(1|3|1)(1|21|1) 


28 


(1|21|1)(2|3|2) 


(1|3|1) (1|21|2) 


(1|3|1)(2|21|1) 


(1|21|1)(1|3|2) 


32 


(1|3|2)(1|21|2) 


(1|3|2) (2|21|1) 


(3)(1|1|1)(1|2|1) 


(3)(1|1|1)(2|2|2) 


36 


(3)(1|1|1)(1|2|2) 


(2)(1|1|1)(1|3|1) 


(1)(1|3|1)(2|2|2) 


(2)(1|1|1)(1|3|2) 


40 


(3)(1|2|1)(2|1|2) 


(3)(1|2|1)(1|1|2) 


(2)(1|3|1)(2|1|2) 


(1)(1|2|1)(1|3|1) 


44 


(1)(1|2|1)(1|3|2) 


(3)(1|1|2)(1|2|2) 


(2)(1|3|1)(1|1|2) 


(1)(1|3|1)(1|2|2) 


48 


(2)(1|1|2)(1|3|2) 


(1|312|1) 


(1|321|1) 


(1|231|1) 


52 


(1|312|2) 


(1|321|2) 


(1|231|2) 


(2|231|1) 


56 


(32)(1|1|1) 


(31)(1|2|1) 


(32)(1|1|2) 


(1)(1|32|1) 


60 


(2)(1|31|1) 


(1)(1|32|2) 


(2)(1|31|2) 


(21)(1|3|1) 


64 


(21)(1|3|2) 


(3)(1|21|1) 


(3)(1|21|2) 


(3) (2|21|1) 


68 


(2) (3)(1|1|1) 


(1)(3)(1|2|1) 


(2) (3)(1|1|2) 


(1)(2)(1|3|1) 


72 


(1)(2)(1|3|2) 


(321) 


(2) (31) 


(3) (21) 


76 


(1) (2) (3) 









Table 1: The basis of (1 2) symmetric six-index polynomials of ki,k2 and rjap- 
Vi = Vi{ki,hi,k2,h2,h3), i = 1, ...,77. We use the short-hand notations: {i\abc\j) = 
ki.ha-hh.hc.kj,..., = ki.kj (abc) = tr{ha-hh.hc),--., (a) = tv{ha) and we omit the 

(1 <-> 2) exchanged term that must be added to all non (1 <-> 2) symmetric polynomials 
of the table. The tensorial dimension of the polynomial Vi is: 5^ = 6 for i = 1, . . . , 15, 
(5t = 4 for i = 16, . . . , 49; 5t = 2 for i = 50, . . . , 73 and (5t = for i = 74, . . . , 77. 
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At criticality we must impose the Ward identities for the Weyl symmetry, which com- 
pare (TTT) traced over one of the stress tensors with the second variation of the anomaly 
( |2.22| ). There are two independent equations by tracing T{ks) or T{ki); in the first case, 
we find (see Appendix B): 

{Q{h){T{ki).hi){T{k2).h2)) 

+ [{{T{-k2)M){T{k2).h2)) + {{T{ki)M){T{-ki).h2))\ 

= A{ki,hr,k2,h2)-\[{h{)A{k2,h2) + {h2)A{kuhi)] . (3.4) 

In this Equation, A{ki,hi) and A{ki, hi; k2, h2) are, respectively, the first and second 
variation of the anomaly ( |2.22| ) w.r.t. to the metric, {hi . 6/6g{ki)), evaluated in flat 
space. Their explicit expressions are: 

A{ki,hi)/X = a' D'^n{ki,hi) = a'kj {kl{hi) - {k^.h^.k^)) ; (3.5) 

and 

A{ki,hi;k2,h2)/\ = (-3c + a) {Rf.ua/s)'^ (h, hi; k2, h2) 
+ (6c -4a) {Tlapf {ki,hi;k2,h2) 
+ (-c + a + r) TZ^{ki,hi,k2,h2) 
+ a' D^n{ki,hi;h2,k2) , (3.6) 

with 

[Rfiuapf {ki,hi;k2,h2) = 2[{ki.k2f'{hi.h2) + {ki.h2-ki){k2-hi-k2) 

-2{ki.k2){ki.h2.hi.k2)], 
{n^pf {kiM-.k2M) = \{kl.k2?{hi){h2) + \klkl{hi.h2) 

+{ki.k2){ki.hi.h2.k2) + (A;2-/i2-^i)(fci-^i-^2) 
+ [\kl{hi){ki.h2-ki) - {ki.k2){hi){k2.h2.kl) 
-kj{k2.h2.hi.k2) + {1 ^ 2}] , 

■R\ki,hi;k2,h2) = 2[{klihi)-iki.hi-ki)){{1^2})], 

D^TZ{ki,hi;k2,h2) = [{kf{hi)-{ki.hi.ki)) 

X {{k1.h2.k1) + {k2.h2.k1) - ^{ki.k2){h2)) + {{1 ^ 2})] 

+{ki + k2)^ [-{hi-h2) (kf + kl + f (fci.fca)) - ^{ki.k2){hi){h2) 
+{k2.hi.h2.ki) + {hi) {{k1.h2.k1) + {k1.h2.k2)) 

+ {h2) {{k2.h1.k2) + {k2.h1.k1))] . 

(3.7) 

The consistency among the two Weyl Ward identities implies the vanishing of the coeffi- 
cient r in the trace anomaly (LI) at criticality; this is the Wess-Zumino condition, which 
requires the commutativity of Weyl variations [l^. Furthermore, the consistencies be- 
tween Weyl and Diff identities imply /q (g^) = /q (A;^) = A a'/3 in agreement with the 
earlier result ( p. 31 ). 



13 



Next, the two Weyl Ward identities are imposed on the form of (TTT) satisfying the dif- 
feomorphism identities, Eq.( |3.3| ); a solution is found in four dimensions only, and the num- 
ber of independent tensor structures reduces from 22 = (17 transverse +5 non-transverse) 
to 13 = (8 transverse-traceless +2 transverse-traceful +3 non-transverse-traceful) . The 
final result for (TTT) off-criticality is thus given by (|3.3| ) up to a convenient change of 
basis, which makes the critical limit simpler and also fulfills other requirements to be 
explained later: 

8 

{{hi.T){h2.T){h;.T)) = Y.Ai{q\k'')% + AE{q\k^)TE + Aw{q\k'')Tw 

i=l 
17 

i=ll 

+ + '^2±] + /a . (3.8) 

± 

In this Equation, the 22 tensors 7^ only depend on the external momenta and do not 
contain any parameter of the theory, thus they are invariant of the renormalization-group 
flow and provide a convenient basis both at and off criticality. The scalar amplitudes Ai 
carry all the dependence on the flow: at criticality, some of them are constrained by the 
Weyl identities, in particular 9 = 22 — 13 of them vanish (7 transverse-traceful +2 non- 
transverse-traceful) . This is made manifest in the basis in Eq. ( |3.8D , which is an extension 
of the basis at criticality. 

More specifically, the 22 tensors have the following properties (see also Table ^): 

• The tensors %, i = 1, . . . , 8 are transverse and traceless w.r.t. all the three stress- 
tensors and thus the corresponding amplitudes are unconstrained at criticality. 

• The tensors Te = Tg and T^y = Tiq are transverse and traceful; their traces match 
the second variation (|3.7D of the Euler and Weyl-square terms in the anomaly ( |1.1| ), 
respectively; the corresponding amplitudes have the critical limits: 

AEiq^k^) - A|(g2,fc2) = A4, 

q^ 

Aw{q'^,k'^) Ay^iq^.k^) = A , (critical points). (3.9) 

q^ 

• The seven tensors %, i = 11, . . . , 17 are transverse and traceful, but their traces do 
not match any of the terms in the anomaly ( |1.1| ). Therefore, their amplitudes should 
vanish at criticality: 

Ai{q^, A;^) ^ , i = 11, . . . , 17, (critical points). (3.10) 
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Table 2: Properties of the tensors in Eq. (|3.8D ; (_L): transversality w.r.t. ki{k2) and 
(yes/no); (tr = 0): tracelessness w.r.t. Ti{T2) and Ta; 5t'- tensorial dimension; A: scale 
dimension of the tensor in polynomial form. 
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These amplitudes are nevertheless different from zero off criticality; another param- 
eterization of (0) off criticaHty has been given in Ref.[Q], but cannot be immediately 
compared with our result. 

The five tensors Ti±, i = 0,2, and 7a are neither transverse nor traceless; the first 
four are the even and odd combinations of the tensors multiplying the two-point 
amplitudes fi (g^) and fi [k"^) (see Eq.(|3^ 



fi± 



1 



{f^ig')±f.{k')) 



0,2 . 



(3.11) 
(3.12) 



The fifth tensor accounts for the one-point amplitude /a ( p.29 ). The critical limits 
of these amplitudes are: 



/0+ {q',k') 
f2± {q\k') 

/a 



A 



-A-^ I los 



(critical points) 



(3.13) 
(3.14) 

(3.15) 
(3.16) 



Actually, the solution of the Weyl identities (3^) yields the results ( p. 13 , 3.14 ), which 
agree with the earlier result ( 2.31 ); on the other hand, the spin- two amplitudes /2± 
are not constrained by the Weyl identities and their critical values are taken from 

(El). 



3.2 Further properties and renormalization 
Null tensor. 

The basis of elementary tensors {Vi} (Table ||) is degenerate in four dimensions by the same 
mechanism found for {AVV) (Eg. (^.41)) and the two-dimensional (TT): there exists one 
"null tensor" which do not manifestly vanish in the tensor notation used here, but actually 
has all components equal to zero in four dimensions. Such six-index Bose-symmetric 
null vector was found by inspection and, moreover, no null vectors were found with five 
and four indices within the tensors bases used for the Diff and Weyl Ward identities, 
respectively. Therefore, the corresponding linear systems for the coefficients have been 
correctly extracted from faithful bases. The null tensor clearly satisfies all Ward identities 
and can be eliminated a-posteriori from the final result (^.§|): it is the traceless-transverse 
tensor Tg (see Table ^), which will be discard hereafter (its expression is given in Appendix 
D). 
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Choice of basis for the tensors. 

Each of the tensors 7^ can be considered as a vector (n) in the 77-dimensionaI Pj-space 
(Table |l]) and the Ward identities are inhomogeneous linear equations of the form A ■ 
u = w, whose solutions possess characteristic ambiguities: (I) the homogeneous solutions 
uq, a ■ uq = 0, are determined up to linear combinations and must be normalized; (II) 
the inhomogeneous solutions are defined up to the addition of homogeneous solutions: 
u = A~^ ■ w A~^ ■ w + Uq. Using this freedom, we can choose a basis of tensors 7^ 
which is more suitable for the physical interpretation. Generically, the tensors have the 
form: 

77 

% = Y. ^^(9'' (^1' ^1' ^2, /i2, /13) , (3.17) 

i=l 

where the Q'f are rational functions of the Lorentz invariants with various scale dimensions 
(the Vi have dimension equal to their tensorial dimension, /S. = 5t = 0,2,4,6, listed in 
Table |l|). Thus, the two ambiguities (I) and (II) can be used to find convenient forms for 
the 

Let us first consider the choice in which the Qf are polynomial: this is clearly possible 
for the homogeneous solutions of the Diff Ward identity (7^, z = 1, . . . , 17), by factoring out 
the common denominators into the corresponding amplitudes. As a matter of fact, four 
of the inhomogeneous solutions (last five entries in Table |^ can also be put in polynomial 
form using type-(II) redefinitions and without changing their amplitudes, which are inputs 
of the Ward identities. The scale dimensions of the 7^'s in polynomial form are listed in 
the last column of Table we find the following cases: 

• The tensors with A = 4 and (5t < 4 are naturally polynomial because they have the 
same dimension as (TTT) ; the corresponding amplitudes are dimensionless functions 
(or constants), as in the case of the two-point function (Section 2.3). 

• The polynomial form also makes sense for the tensors with A = 6 and 5t = 4, 6; 
their amplitudes have the natural dimension —2 and may display a simple pole in 
the massless theory. 

• The polynomial form is not convenient for tensors with (Jj- = 6 and A > 6; the 
corresponding amplitudes would have large negative dimensions (4 — A = —4, —6) 
which suggest the possibility of unphysical higher-order poles to be cancelled by 
corresponding zeroes in the tensors. By inspection, we have found another non- 
polynomial form for these tensors which is characterized by scale dimension 6 and 
rational Qf with simple poles ~ l/('7^ — 2/c^), l/{q'^ + 2/c^) (higher-order poles 
and l/k"^ poles are absent in this basis). The corresponding amplitudes have again 
dimension —2. 
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In conclusion, the general (TTT) can be put into a form which does not explicit contain 
higher-order poles and possess a smooth limitp] k"^ ^ with q"^ 0. This form is convenient 
in the kinematical region > 2k'^ > 0; other choices of basis may be relevant for different 
corners of the phase space, say q^ ^ with k'^ / 0. 



Renormalization and scheme-dependence of the amplitudes. 

As discussed at the end of Section 2.3, the form of (TTT) in ( |3.8| ) can be changed by 
adding counterterms, which are polynomials of the momenta of maximal dimension four. 
The counterterms cannot be arbitrary, because they would change the definition of the 
stress tensor, e.g. by transforming a trace anomaly into a gravitational anomaly; therefore, 
the counterterms should be already present in the solution of the Ward identities just 
found. We can identify them as the tensors with A = 4 and 5t < 4, which actually are 
polynomials; the corresponding dimensionless amplitudes are scheme-dependent by the 
shifts — > + const.. Such tensors are: 

• The null tensor Tg which vanishes identically. 

• The tensor 7o+ which matches the scheme-dependent D^TZ term in the trace anomaly; 
this verifies the ambiguity a' — > a' + const, at the level of the three-point function. 

• The tensor 72+, which accounts for the /^-dependence of the logarithms in the ampli- 
tude /2+ (Eq.( |3.15 )); such dependence instead cancels out in the odd combination 



/2_, which multiplies a A = 5t = 6 tensor. 
• The tensor T\ coming from the one-point function off criticality. 

These results are in agreement with earlier analyses based on the lower-point functions 
and the effective action and provide a check of our derivation. 

Note that two further tensors, Tn and 7i2, have dimensions = 4 and A = 6, so that 
they cannot mix with counterterms. Here we have checked that their dimension cannot 
be further lowered by changing the basis: the mixing with other tensors is not actually 
possible owing to their higher 6t value. 

In conclusion, the 6t = Q tensors in Table |2| (and the previous Tn and T12) identify 
amplitudes Ai{q'^, k'^) which are finite and scheme-independent with scale dimension (—2). 
In particular, we found the desired result of matching the (universal) anomaly coefficients 
a, c with two of these amplitudes, namely Ae and Aw in (|3.9D. 



^ Note that all tensors become polynomial at = 0; the solution of the Ward identities at fc^ = will 
be further discussed in Section 4. 
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5 
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1 
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10 
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1 
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8 
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(5t = 4 
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Table 3: Summary of the properties of the solutions to the Ward identities for k"^ and 
k"^ = 0. Using the notations of Table ^, we report the number of solutions, split by the 5t 
value, both off criticality (Diff line) and at criticality {Weyl); we also indicate the critical 
limit of the amplitudes (value) and their labels. 



Limit of the solution for k"^ ^ 

As already said, the limit is smooth in our choice of 7^ basis; however, one tensor becomes 
a linear combination of the others and should be removed. Actually, for k"^ = 0, the 
two tensor structures of the two-point function {{hi.T{k)) {h2.T{—k))) are degenerate 
(see Eq.(^30|)), so that the Diff Ward identity looses one inhomogeneous solution. 



Moreover, the amplitude /2 [k"^) becomes singular at criticality (see Eq.( |2.32|) ); thus, the 
corresponding term /2 (A;^) [72+ — ^2-] is subtracted from (TTT): 

{{h,.T){h2.T){hs.T)) 

reg, k-^=0 

= hm^|((/ii.r)(/i2.r)(/i3.r))-^^^[r2+-T2_]| . (3.18) 

The k'^ = solution thus contain 21 = (17 transverse -|-4 non-transverse) independent 
tensors. The analysis of the Ward identities discussed before has been independently 
repeated at the point k^ = and matched with the result of the limit from k'^ ^ 0. A 
summary of the properties of the solutions for both A;^ = and fc^ 7^ is given in Table 
^; note that some of the tensors have lower St values at fc^ = 0, but the number and type 
of homogeneous solutions remains the same. 
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3.3 Critical points: conformal isometries of flat space 



In the previous Section, we obtained the general form ( |3.8| ) of (TTT), which contains 9 
scheme-independent ampUtudes at criticahty, namely Ai,i = 1, . . . ,7, Ae and {A]y, /2-). 
Previous approaches in coordinate space have instead found three independent non-contact 



tensor structures |19|[^0[: two parameterized by the a and c anomalies and one traceless 
transverse". The disagreement is due to the fact that, in our derivation, we did not 
impose the symmetry under the regular (i.e. global) conformal transformations of flat 
space, belonging to the group SO{d + 1,1) |l^|jl5|; as is well known, these symmetries 
almost completely determine the functional form of the two- and three-point functions 
in coordinate space. Moreover, we did not considered the functional relations among 
amplitudes due to the exchange symmetry of T(ki) with T{k^). 

On the other hand, the disadvantage of the coordinate-space approach is that it yields 
the expressions of correlators for non-coincident points, which are rather difficult to renor- 
malize Jl9| . Furthermore, the coordinate approach cannot be easily extended away from 
criticality because the Diff Ward identities become involved. 

In the following, we rederive the Ward identities for the SO{d +1,1) conformal sym- 
metry in flat momentum space, without neglecting the contact terms and consistently 
following the conventions set in this paper. As is well known jl^, a SO{d + 1,1) con- 
formal transformation is a combinations of a reparametrization and a Weyl scaling which 
map the Minkowski metric into itself: the precise notion is that of conformal isometries 
of the flat metric, which are associated to conformal Killing vectors [23]. A general metric 
g^^{x) possesses a conformal isometry if it satisfies: 

= (5e + 5„) g^"'{x) = D(^'e^)(x) - 2a{x) (7^"^(x) , (3.19) 

where is the conformal Killing vector; the associated Weyl scaling cr{x) = {D ■ e)/d 
is found by taking the trace of ( |3.19| ). In flat space, the solutions of ( |3.19| ) generate the 
SO{d + 1,1) conformal group |jl^ 



e^'(x) = b^' + Xxf" + Loiix" + a^x^ - 2(a • x)x^' , {d > 2). (3.20) 

We recall that = —^e is (minus) a Lie derivative and that the combination (6^ + 6^) 
is still covariant because the Weyl transformation is local. The usual conformal Ward 



identities in flat space |15] are nothing else than the covariance conditions for the flelds. 



For example, a scalar and a vector fleld satisfy: 



" The three independent structures were originaUy found in free field theory by the first of Refs.|19| 
and were confirmed by the subsequent analyses of the conformal symmetry. 
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A- 1 

d 



{d-e)A, 



(3.21) 



The field (p, A^ are said quasi-primary, because they transform homogeneously under Weyl 
with scaling dimension A. In the case of the stress tensor, we must pay attention to the 
effects of the trace anomaly, and the naive (classical) formulas ( 3.21| ) should not be taken 
for granted. 



We reconsider the Ward identities ( 2.20 ) for Diff and Weyl transformations of the 
generating functional and specialize them for a conformal isometry ( 3.ig| ) (for more details, 
see Appendix B). The two identities can be combined to obtain the sum of three terms: 
respectively, the stress-tensor (?7- + l)-function times {5^ + 5a)g, the n-th point function and 
the n-th variation of the anomaly. Since the first term vanishes for conformal isometrics, 
this Ward identity closes on the n-point functions and yields a finite set of differential 
equations with inhomogeneous term given by the anomaly. These equations put additional 
constrains on the solutions of the two separate Ward identities, which apply to general 
metrics without conformal isometrics. 

The conformal Ward identity for the two-point function is: 



[(5, + 5c.)i + {5e + 5a)2] {Taa'ixi)Ti3i3>{x2)) 



-2 



5g'^»' {xi)6gl'f^' {X2) 



Vg o- A{g) 



(3.22) 



a={d-e)/4 

where the metric variations are evaluated in flat space. The conformal transformation of 



the stress tensor reads |15| : 



Tap + [d^^ae") T,, 



/3) 



(3.23) 



For scale variations (e^ = b x^^) and in four-dimensional momentum space, Eq.( 3.22| ) 
becomes: 

{p-dp- 4) {{hi.T{p)){h2.T{-p))) = 2 A{p, hi; -p, h2) , (3.24) 

where the r.h.s. is the second variation of the integrated trace anomaly ( |1 . 1] ) in the notation 
introduced before (see Eq. (|3.7| )). This variation only contains the Weyl term, because the 
Euler term is topological invariant and ^D^TZ is a total derivative; therefore, the r.h.s. 



of (3.24) is a traceless-transverse tensor. By plugging in the explicit expressions of the 
anomaly (|3.7f ) and of (TT) (|2.30| ), one verifies that the coefficient of the logarithm in 
spin- two amplitude /2 (p^) is equal to the anomaly number c, as anticipated in Eq.( 2.32| ) 



For special conformal transformations, parameterized by in (3.20), the Ward iden- 
tity (|]2|) becomes a second-order differential equation in momentum space; we do not 
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solve it but rely on earlier studies in coordinate space (see e.g. Ref. ||5|) to conclude that 
it is automatically satisfied by the critical form of (TT) given by ( |2.3C| ) with /g , /2 in 
(|3l|,|]3|). 

In conclusion, we have found that the conformal Ward identities for the stress-tensor 
n-point correlators contain an inhomogeneous term given by the n-th variation of the trace 
anomaly, and that they amount to a set of differential equations for the scalar amplitudes 
at criticality. At the three-point level, the equations for special conformal transformations, 
together with the functional relations for the (1 3) symmetry, presumably determine 
the form of the traceless transverse amplitudes, Af{q'^,k'^), i = 1,...,7, leading to an 
agreement with the coordinate-space result [|19[||2[1U. We shall not carry out this analysis, 
because we are more interested in off-critical behavior of the traceful terms in (TTT), 
which describe the renormalization-group flow of the anomaly coefficients a and c. 



4 Trace anomaly: reduction formulae and sum rules 
4.1 Reduction formulae 



We now describe the inversion of the (TTT) expansion ( |3.g| ) which extracts the amplitudes 
^a;(Q'^7^^) from an explicit expression of the three-point function, such as the result of a 
calculation. The latter can be expanded in the canonical basis of Bose-symmetric tensors 
Vi, i = 1, . . . ,77 (see Table |^), leading to a set of know coefficient functions Ui{q^, fc^). We 
then write (see also ( 3.17D ): 
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(TTT) = ^ A^T^ = ^^ ^..((/^A;2) Q-(g2,A;2)P,(A:i,/ii,A;2,/i2,/i3) 

U)=l 1 = 1 L0 = 1 

77 

= uM\k'')V^. (4.1) 

1=1 

After the elimination of the null vector, the 21 7^'s and the 76 ViS form two sets of 



linearly- independent "vectors"; therefore, Equation (4.1) defines the linear problem of 



finding the "coefficients" by inverting the 21x76 matrix Qf of maximal rank. The 
solution is unique and expresses the 21 j4^'s as linear combinations of the Uj's, together 
with 55 linear relations among the UiS. 

These constraints verify that the calculation, besides being correct, satisfies the Diff 
Ward identities in our specific formulation, which includes a complete choice of renormal- 
ization scheme and of stress-tensor definitions. This leads to the practical problem that 
a given calculation may not immediately fit our scheme, and thus the check of our Ward 
identities (i.e. of the 55 constraints) may require modifications of the calculation rules and 
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become a painful process. For example, in Section 5 it is shown that a Feynman-graph 
calculation of (TTT) needs the inclusion of T^^/ vertices with an increasing number of legs 
in order to match our definition of the stress tensor by differentiation of the generating 
functional, Eq.( ^.19| ). There is basically no short cut to this problem if one really needs 
to determine the whole expression of {TTT) including the scheme-dependent terms. 

An important simplification occurs in deriving the scheme-independent 5t = Q ampli- 
tudes: firstly, one can limit the inversion ( [4.1| ) to the 15-dimensional subset of 5t = 6 
basic vectors Vi, i = 1, . . . , 15 (see Table ||); secondly, these amplitudes are independent 
of renormalization scheme and T^,y definition^, so that eventual constraints among the 
tij's, should be automatically verified in any calculation[^. More precisely, for k"^ the 
number of amplitudes is 14-1-2 (transverse-|-non-transverse) (see Table the two non- 
transverse tensors are not linearly independent in the 5t = Q basis but otherwise known, 
thus there are 15 unknowns and no constraints. For k"^ = 0, there are 12-1-2 unknown 
amplitudes (the inhomogeneous ones are still dependent) and 2 constraints; finally, at the 
critical point there are 9-1-1 unknowns and 5 constraints in both cases. 

Among the 5t = Q amplitudes, we are particularly interested to Ae and A\y, which 
reduce to the anomaly coefficients a, c at criticality (see Section 3.1). Hereafter, we present 
the determination of the St = Q amplitudes at the point /c^ = 0, which is sufficient for our 
later purposes (the general case can be found in the Mathematica notebook associated 
to this paper). The 12-1-1 independent A;^ = amplitudes are obtained by inverting the 
(numeric) 13x15 matrix Qf{q^,k'^ = 0) of Eq.( |4.1D restricted to the subspace of {6t = 6) 
Pj's. We remark that this linear problem can be equivalently solved by introducing a dual 
basis of tensors T*, {T*,T^/) = 5^^^}', where the brackets denote a suitable degenerate 
scalar product which vanishes on the 5t < 6 basis Q. 

We finally obtain the following pair of reduction formulae for the anomaly amplitudes 
(A:2 = 0): 

^E(g^O) = ^{4u2-4u4-U5 + 2u8-2uu + 12uu + 3ui5) , (4.2) 
3d 

Aw{q^,0) = (-Sua + 8^4 + 2^5 - 4?/8 + 4^11 + 12^14 + 3ui5) , (4.3) 

lUo 

together with the two constraints: 

Ul = ^U3 , ui2=0 ■ (4.4) 



Note that the formulae ( |4^j4^ ) are invariant under the changes of basis, T^ — > 7^ -|- 

Yli=i "^i^i^ ^ = E, W, which add traceless-transverse terms to the anomaly tensors. 

** An argument for the definition-independence is given in Section 5. 

Note that the null vector has 5t = 4 and is automatically excluded from this basis. 
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Another interesting reduction formula yields a combination of the (TT) amplitudes 
which is associated to a, 6t = Q (TTT) tensor (see Table ^): 

^ (2/2(9^) + fo{q') - Uk' = 0)) = ^ ; (4.5) 

this result gives a useful check. The complete set of 13 reduction formulae is reported in 
Appendix D. 

At the critical point, the number of independent tensor reduces to 10 and there are 5 
constraints; in this case, the reduction formulae determine the anomaly coefficients (see 
Eq.(IJ)): 

a = (-""11 - ^13 + 4wi4 + ni5) , (4.6) 

c = (2uii + 2ni3 + 4^14 + U15) ■ (4.7) 

The constraints at criticality are: 

■ui = ^ M3 , -"12 = , 

U2 = ^ (4n4 + Us - 2us - un - 8^13) , 

""3 = ^ (9^6 - 2^11 - 2ui3 - 4^14 - ^il5) , 

y 

tig = ^ (-uio + nil + ni3) • (4.8) 

Equations ( [4.6| , p7| ) allow to compute the trace anomalies by standard flat-space pertur- 
bative calculations, and can be useful for interacting theories. 

Infrared divergences. As emphasized, the amplitudes Ui{q'^,k'^ = 0),i = 1, . . . , 15 are 
ultraviolet finite and satisfy un-subtracted dispersion relations of the type ( |2.14| )) away 
from critical points. However, some amplitudes may develop infrared singularities at 
criticality, which are introduced by the kinematical limit A;^ = 0. This fact is illustrated 
by the formula (|4.5| ) (see also Table ^). For k'^ 7^ 0, the {6t = 6) tensor T2- would identify 
a finite amplitude with safe critical limit m ^ 0, namely /|„/(7^ ~ log{q^ / k'^) / . For 
k'^ = 0, the critical limit is instead /K^^)/^^ ~ \og[q^ / m?) / q^ , with the physical mass 
m acting as infrared regulator (rather then the renormalization scale /i controlling the 
ultra-violet singularity of the (TT) amplitude ( 2.32| )). 

The standard solution to this problem amounts to the choice of another point A;^ 7^ 
for the reduction formulae at criticality. However, this is not necessary for the derivation 
of the anomaly amplitudes Eqs.(4.2,4.3,[4.6,4?^), which are free of infrared singularities. 
Actually, these amplitudes are related to the trace of the stress tensor, which vanishes 
as an operator at criticality due to Weyl invariance; thus, these amplitudes contain the 
anomalous pole only. 
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4.2 Sum rules for the renormalization-group flow of a and c 



We start by recalling the sum rule for the two-dimensional central charge c, Eq.( |1.2D [18|: 
we present a simplified version of the dispersive proof of Ref . ^] , which can be conveniently 
generalized to four dimension. In Section 2.3, we wrote the general form of (TT) satisfying 



the DifF Ward identity, Eg. ( 2. 30 ); in two dimensions, this reduces to: 

iih.Tip)) (h2.T{-p))) = f (/) {p\h^) - {p.h-p)) {p\h2) - ip.h2.p)) . (4.9) 

The spin-zero scalar amplitude f{p^) has dimension (—2) and is finite and scheme- independent, 
namely its renormalization is analogous to that of the four-dimensional amplitudes vis ((7^, k'^) 
of (AVV) (Section 2.1) and Ae,A\y of (TTT); note, however, the difference in the single- 
variable dependence. The critical limit is: 

/(P^) ~^ o ) (critical points), (4-10) 

3 

where c is the Virasoro central charge and T has the standard two-dimensional normal- 
ization Q. The amplitude satisfies the un-subtracted dispersion relation: 

_ f * W(£) _ „ jl, 

J IT S — p^ 

Im f{s) = ^ c{s) , (4.12) 



which also identifies Im /(s) with the spectral measure c{s) introduced in ( |1.2[) . Its critical 
limit is: 

c{s) c 5{s) , (critical points). (4-13) 

We now introduce the following sum rule off-criticality, 

Y.{m^) = — dslm f(s,m^) = ds c(s,m^) , (4.14) 
Jo Jo 

and explicitly write the dependence on the mass scale m introduced by the relevant (renor- 
malized) coupling constant g. In the following, we would like to prove that $](m^) is 
actually independent of m, 

S(m2) = S(0) , for < < 00 . (4.15) 



From Eg. ( [4. 13 ), we know that S(0) = cjyy, the central charge of the ultra-violet fixed 



point. If S is independent of m, it is also renormalization-group invariant, because t ~ 
log(m) is the parameter of the renormalization-group fiow between the fixed points (UV) 
and (IR). 
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Equation ( 4.15| ) can be derived from the general property that the (trace) anomaly is 



at most a constant; this was shown in the case of {AVV) by using dispersion relations 
(see Eqs.( p^H2l7|) ): since the imaginary parts of amplitudes satisfy non-anomalous Weyl 



Ward identities, the anomaly should be accounted for by the subtraction constants. How- 
ever, the precise formulation of this argument for (TT) (and later on for (TTT)) requires 
the dispersive analysis of all the amplitudes related by the Ward identity; we prefer to 
use another argument based on the dimensional analysis and the renormalization group, 
which is presented hereafter. 



The amplitude fijp ) in (|4.g| ) is finite and scheme- independent: there is no wave- 
function renormalization because the stress-tensor has no anomalous dimension, while 
the coupling-constant renormalization makes the quantity finite and brings in the scale 
m. Infrared divergences are also excluded for this quantity in two dimensions Q. The 
imaginary part of f{p'^,m?) has support on the real positive s-axis and goes to infinity 
as o(l/s), i.e. faster than the behavior allowed by its scale dimension, such that the 
integral in ( [4.14] ) is finite. This extra convergence can be verified perturbatively close to 
the UV fixed point since /(p^,m^) measures the trace anomaly, it corresponds to 

an expectation value involving one O operator (at least); we thus can write, 

f{p\m^) = {@{p)...) = [3{g) im---) 

~ y g s^^^"^^^ , for s — > oo , (4.16) 

where y > is the dimension of the coupling constant (go = g m^), and we used the 
relation Q = f3{g) (p expressing the trace off-criticality as the beta function times the 
perturbing field [p^. 

In conclusion, S(m^) ( [4.14| ) does not contain any regulator for the integrand or the 
integral; it depends on a single dimensionful parameter, i.e. m, and thus is necessarily 
constant (and renormalization-group invariant). 



The sum rule (4.14) implies some conditions on its integrand c(s); we can write, in 
general Q: 

c(s) = Co 5(s) + ci(s,m^) . (4.17) 

The function ci{s,m?) should be a representation of the delta- function with spreading 
parameterized by m^, so that the ultra-violet limit m — > matches the critical form 
( 4.13|) . In the opposite limit m —>■ oo, ci(s,m^) becomes very broad and it contributes 



infinitesimally to the correlator (TT) for any fixed scale s = sq; therefore, c{s) must be 
redefined in the infrared scale- invariant theory (m = oo) by means of the weak limit: 

c(s) ^ Co (5(s) , Co = c/R , (m = oo). (4.18) 
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This Equation identifies cq with the infra-red value of the central charge (cq is another 
renormahzation-group invariant). We can finally rewrite the sum rule ( [4.14] ) in the more 
famihar form ([lD [|18|: 



CUV 



/ ds c{s) + I ds c{s) = ciR + / ds ci{s,m?) . (4.19) 

JQ Je Je 



This equation implies the c-theorem in the integrated form, cjjv > c//j, because the 
spectral measure c(s) is positive definite ||l^. Furthermore, the integral of c(s) up to a 
fixed (finite) scale defines a function c{g) off-criticality |5|, which interpolates between 
cijv and cjji and satisfies the c-theorem in the differential form [^i dc{g) / d log (m) = 
-(3{g) dc{g)/dg < 0. 

The four dimensional sum rules for the anomaly coefficients a and c are obtained by 
similar arguments. The corresponding scalar amplitudes A^{q'^,k'^ = 0), lo = E,W, are 
extracted from the three-point function (TTT) by the procedure explained in the previous 
Section; this also applies to the imaginary parts w.r.t. because the tensors %j are 
polynomial at A;^ = 0: 

Im {TTT) = Im As(s, 0) Te + Im Aw{s, 0) Tw + ■ ■ ■ , {k^ = 0). (4.20) 

The amplitudes have dimension (—2) and satisfy un-subtracted dispersion relations; their 
critical limits ( p.9| ) imply: 

Im Ae{s,Q) — vrA a 6{s) , 

Im A\Yis,0) —7rXc6{s), (critical points). (4-21) 



The sum rule for the a anomaly is defined in analogy with ( 4.14 ): 



1 /"°° 

T.E (m^ O) = / ds Im Ae (s, k'^ = 0, m^) , (4.22) 

■^■^ Jo 

with the m-dependence explicitly written. This integral is meant to extend over the 
support of Im Ae- 



We now follow the same steps as for the two-dimensional case, Eqs.(|l|-|l|): we first 



argue that S_e is m-independent and equal to the ultra-violet value ajyy, provided that it 
is finite and does not involve any regulator. The integrand is indeed finite and scheme- 
dependent, as already stressed; the integral is convergent if the stress tensor is "improved" 
at both the UV and IR fixed points, i.e. = at criticality (apart from anomalies). Ac- 
tually, Im Ae is related to an expectation value involving the G operator, whose vanishing 
at the UV (IR) critical point forces the asymptotic behavior of Im Ae{s,0, m?) to be faster 
than 1/s, the scale-invariant law, for s ^ oo (respectively, s — > 0). 
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The improvement is implied by Weyl invariance, assumed throughout this paper, but is 
not guaranteed for a generic fixed point: as already said, four-dimensional scale invariance 
only require G to be a derivative operator. Near a generic UV fixed point, it is rather 
natural to assume the improvement and the relation G = f3{g) (j), Eq.( [4.16 ) [28|. Near the 



IR point, it is less obvious: for example, the spontaneous symmetry breaking of a global 
symmetry may lead to an infrared critical theory without improved stress tensor We 
are aware that the hypothesis of improvement at the IR fixed point (convergence of the 
sum rule) may exclude interesting non-perturbative renormalization-group flows from our 
analysis. 



The constancy of the sum rule ( 4.22| ) again constrains the amplitude to be of the form 
&■ 

-J^IuiAe (s, k"^ = 0, m^) = aMs) + ai(s, v?) ; (4.23) 
ttA 

where the smooth function ai is a representation of the delta function for m ^ 0, and 
the explicit flow-invariant delta term survives in the infrared theory. Note that the scale- 
invariant power law, Im A£;(s, 0, 0)(is ~ ds/s, is excluded at criticality because it would 
imply a non-vanishing G operator and a divergent sum rule, as emphasized. 

We finally obtain the sum rule for the four-dimensional anomaly, analog of Eq. (|1.2|): 

1 f°° 

ajjv — o-iR = / ds Im AE(s,0,'m'^) . (4.24) 

ttA J, 

The other four-dimensional anomaly number c satisfies a similar sum rule involving Im A\y{s, 0, m? 
Equation ( [4.24| ) is the announced result of this Section, which is followed by a number of 
remarks: 



If Im AE{s,0,m'^) were positive definite, then Eq.( [4.24j) would prove the c-theorem 



in four dimensions, auv > CLiRi ^ suggested by the non-trivial examples mentioned 
in the Introduction [^][^. We do not discuss the positivity in this paper, but only 
add some remarks in the Conclusions. 

On the other hand, the corresponding result for the Weyl-squared term, cjjv > cjr, 
is not true in general, as shown by known counter-examples [Q|6|; the respective 
amplitude Im A\y{s, 0, m?) cannot be positive definite in general, although its critical 
limit is always positive, i.e. cuviCm > 0, due to the positivity of (T(x)T(O)) at 
criticality (|x| / 0) |p. These facts suggest the absence of naive back-of-the-envelope 
positivity arguments. 

The other three 6t = Q amplitudes Ai corresponding to traceful tensors % {i = 
15, . . . , 17 in Table y) also satisfy analogous sum rules; however, their vanishing 
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critical limit imply null sum rules ( 4.22 ), Sj = 0, for i = 15, . . . , 17: thus, the 
imaginary parts of these amplitudes cannot have definite sign. 



We stress that the positivity of Im Ae{s, 0, m?) is not necessary for the derivation of 
the sum rule ( [4.24 ); its smooth part ai{s, m?) can be any representation of the delta 



function which becomes very broad in the infrared limit m ^ oo (this is assured by 
the decoupling of massive states). 

We conclude this section with a remark on the chiral anomaly, analysed in Section 2.1: 



the sum rule ( 4.24| ) can also be written for the flow of the chiral anomaly, as is clear from 



Eg. ( 2.1?] ), analogous to ( 4.22 ). However, there is an important difference: the role of the 



Weyl symmetry is now played by the chiral symmetry, such that "critical theory" should 
be replaced with "theory with massless fermions". These are not equivalent properties, 
because it is conceivable a renormalization-group flow with chiral symmetry, for example 
a QCD-like theory with massless quarks. In this case, the relevant operator causing the 
flow yields the stress-tensor trace according to ( |4.16| ), but does not couple to the chiral 
current; as a consequence, Im ^3(5,0) in ( |2.17| ) remains equal to its UV critical form of a 
delta function and there is no actual flow of the chiral anomaly [ p^ : 

Cuv = CiR , (RG flows in chiral theories). (4-25) 

Furthermore, the term Im {AVV) ~ 5{s) means that the chiral anomaly is saturated by 
massless on-shell intermediate states {q^ = /c^ = 0) ||3^; thus, it can be perturbatively 
computed either from the fundamental (UV) or the composite (IR) massless fermions; 



Equation ( 4.25 ) corresponds to the 't Hooft anomaly matching conditions |£2[. In conclu- 
sion, the sum rules for the conformal and chiral anomalies are mathematically identical 
but may lead to rather different physical consequences. 



5 Test calculations in free field theories 

We now apply the results of Section 4 to the cases of the free massive scalar and Dirac 
fermion theories[]; we (re)compute here the anomaly coefficients a, c, verify the related 
sum-rules and check other properties of the reduction formulae. 

The stress-tensors three-point correlator is obtained, according Eq.( |2.19| ), by triple 
variation w.r.t g^^ of the generating functional: 

W{g) = log j [d(t)\e^^^-^3) . (5.1) 
^Background material for these calculations can be found e.g. in Ref. H]. 
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After variations and flat limit one schematically obtains: 



(TiTaTg) = {5il 82! 62l + {8il 826:^1 + perms.) + 5i525zl) + ---, (5.2) 

where we denoted by perms, the permutations of labels and we left unspecified some 
disconnected terms, all containing a {51) factor. In terms of Feynman diagrams, the 
first three contributions to (TTT) in Eq. (5.2) are described respectively by three one 
graviton vertices simply linked in a triangle, by a two graviton vertex doubly linked to a 
one graviton vertex, and a tadpole over a three graviton vertex. 

It is clear that a complete calculation of all contributions - even in a free field theory - is 
a nontrivial algebraic task. Nevertheless, it should be reminded that the results of Section 
4 only involve the subset of (TTT) terms that have tensorial dimension 5t = 6. We now 
remark the following simple property: if the action / contains at most two derivatives, any 
insertions of (variations of) / in the expectation value can raise the tensorial dimension 
by two at most. It then follows that the second and third contributions to {TTT) have 
tensorial dimension Jj^ < 6 and should not be considered for our purposes. Moreover, 
owing to the fact that {Sgt^^I) oc d/^n, in a flat space, the leftover disconnected terms in 
(|5.2[) also have 6t < Q- 

In conclusion, we are left with the triangle diagram whose explicit contribution is, 

/d'^L 
j^Vh,{L - q, -L)S{L)Vh,{L, h - L)S{L - h) 

xVh,{L-h,q-L)S{L-q) , (5.3) 

where the vertex V{ki,ki) and the propagator S{p) can be found, e.g. in Ref.|^](see 
Equations (3.27-3.30) there). 

The imaginary part w.r.t. of (rrT)triangie for q'^ > and /c^ = is obtained by 
the standard Cutkowsky rules, cutting the lines linked to the vertex T3 (the remaining 
integrals are elementary). After collecting from the resulting expression the coefficients 
Ui{q'^,0) of the expansion in the canonical {6t = 6) basis Vi, i = 1, ... 15, we have tested 
the formulas of Section 4. We have found that: 



The constraints (4^) and (|4.8|) are satisfied off and at criticality, respectively. 



The reduction formulae (4^,4^) determine the following expressions for the Euler 



and Weyl anomaly amplitudes, in the massive scalar theory: 



— Im Ae{s, 0,171^] 

ATT 

e(x - 1) — 



32rx(l + 3x) + (-27 + 80x + Ux'^) log (- — -] 

\1 + 



(5.4) 
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— Im Awis, 0, m^) 
Avr 

-9{x-l) 



8rx(13 + 3x) + (27 + 26x + 12x^) log 



1 -r 
1 + r 



with 



Am? 



In the massive fermion theory, we find: 



- — Im Ae{s, 0, m^) 

ATT 



24^2^,4 



8rx(4 + 17x) + 3(-9 + 25x + 12x^) log 



1 - r 
1 + r 



1 



Im Aw{s^ 0, 
5 



Avr 

efx-l 



72m^a;^ 



4rx(52 + 5x) + 3(18 + 23x - 2,x^) log 



1 — r 
1 + r 



(5.5) 



(5.6) 



(5.7) 



(5.^ 



The sum rules (4.24) for these imaginary amplitudes verify the expected renormalization- 
group fiows of the trace anomalies (see e.g. Ref. for the well-known free-field 
anomalies): 



auv — aiR = 1 , CUV — ciR = 1 , 
auv - aiR = 11 , CUV - ciR = 6 



(scalar); 
(Dirac fermion). (5.9) 



The amplitudes Ai, i = 15, 16, 17, of the traceful tensors vanishing at criticality have 
also been computed (using the formulas in Appendix D); the corresponding sum 
rules ( 4.22| ) are convergent and equal to zero, as expected: 



S,: = 



15,16,17. 



(5.10) 



The imaginary parts of the Euler and Weyl amplitudes turn out to be positive 
definite in all the above examples but one, i.e. the bosonic Im Ae ( |5.4| ); nevertheless, 
a positive definite function can be obtained by adding to Im Ae a suitable linear 
combination of the three non-definite amplitudes Im Ai, i = 15, 16, 17, which have 
vanishing sum rule. 

Infrared divergences occur at criticality for some of the amplitudes Ai, i = 1, . . . , 7, of 
the traceless tensors, in agreement with the discussion in Section 4.1; moreover, the 
formula f2{<l^)/Q^ ~ \og{q^ / m?) / is checked for m ^ (including the numerical 
coefficient computing c again). 
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We finally remark that another possible test of the reduction formulae at criticality 
( ^■6| , [4.7]j4^ ) is provided by the free massless vector field; however, the theory becomes 
interacting off-criticality (there is no flow into the massive free vector), so that the sum 
rules can only be checked perturbatively. We postpone this analysis to further develop- 
ments of our work. 



6 Conclusions 

In this paper, we have found the general four-dimensional form of the stress-tensor three- 
point function, written as an expansion over scalar amplitudes times tensor structures 
solution of the Ward identities. We thus made explicit the conditions imposed by the 
mandatory symmetries, leaving out any dynamical hypothesis. From this result, we have 
identified the amplitudes corresponding to the two universal trace-anomaly numbers a and 
c, and derived the renormalization-group sum rules which they satisfy. 

The main motivation for this analysis has been the extension of the Zamolodchikov c- 
theorem to four dimensions: we did not provide a proof but we tried to clear the ground for 
further attempts to it. One possible road could be the proof of positivity for the imaginary 
amplitude entering the a sum rule ( ^.24| ) (possibly using the weak-energy condition of 



Ref.[|0|). 



Other arguments which stands on specific dynamical hypotheses can hopefully take 
advantage of our reduced set of variables. For example, exploiting the conditions of = 1 
supersymmetry, which is present in most of the non-trivial tests of the would-be a-theorem 



It would also be interesting to reformulate the dynamical assumptions of Refs.|LO| 
at the level of three-point functions. Finally, the analysis of manifestly positive-definite 
four-point amplitudes is feasible with the technology presented in this paper. 
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A Notations and conventions 

The four-dimensional trace anomaly is written: 

9^''{T^u)^{e) = X{aE-3cW + a' D^n + rTZ^) , A = ^^^^^ . (Al) 
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where the Euler density and the square of the Weyl tensor are, respectively, 

and is the square of the covariant derivative. The Riemann and Ricci tensors are: 

^2^3/^4 ^ ^M3-^/i4M2 ^t^i ^J■3^J■2 ~ IJ-3IJ-5 l^4fJ-2 ^4/^5 /^3^2 ' '^^-'^^ 

^M2/.4 = ^;:w4' (a.4) 

The LeviCivita connection is: 

''^/^2At3 ~ 2^'^^^^ i^^J■2 9^lA^^z +^M35At4A'2 ~ ^M4S'At2/i3) ' (^-5) 

and the metric has Euchdean signature. We also use the Minkowski signature (+, 

— ) when discussing dispersion relations (Sections 2.1 and 4.2); on this metric, the trace 
anomaly ( [A.l| ) changes sign and the normalization A should be chosen negative. 

For any local quantity A{x]g) functional of g'^^'{x), we consider its n-th variation at 
the point and contract the pairs of indices with generic polarization tensors 

hf^ as explained at the beginning of Section 2.3. We denote the n-th variation by: 

A{xuhv,...;xn,hn;x)^ h't'^'- J? • • • ^ A{x- g) . (A.6) 

The differential of a functional F[g] w.r.t. the metric, g'^^' g^^' + A^^', with A an 
infinitesimal function, is: 

<5a F[g] = F[g + A] - F[g] = J dxj-^^ A^^'(x) . (A.7) 

The Fourier transform of translation-invariant functions f{xi,X2, .., a;„) is defined by: 

f{xi,X2,..,Xn) = (27r)'^""'^y" dpi...dpn 6 (Y^P^ (^Y^P"-'^^ -^(^1' • 

(A.8) 

In particular, the metric variation ( [A.6| ) is translation invariant after flat-space limit and 
its Fourier transform is denoted hy A {ki, hi; kn, hn) ■ 



B DifF and Weyl Ward identities 

The Diff Ward identities for the stress-tensor n-point functions are obtained by performing 
the reparametrization 5^ on the argument of the generating functional 14/^[gf + /i]; here, g^'^ 
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is the background metric and h^^'^ is a small deformation, which is used as a source for 
extracting the stress-tensor correlators according to (2.19). The Diff transformations are: 



(B.l) 
(B.2) 



Upon expanding 1^[(7 + /i] in series of h and performing the transformations, one finds 
that each term 0{h'^), n = 0, 1, . . . yields the Ward identity for the (n + l)-point function 
(also involving the n-point correlator); the local form of the identity is obtained after 
integration by parts of the derivatives in ( B.1 , |B.2| ) and the limit to flat space. To order 
0{h^), one finds d^{T^i,) = 0, and to order 0{h) the two-point identity Eg. ( p. 23 ). The 



three-point identity (3.2) becomes simpler in momentum space and after contraction of 
indices with polarization tensors, as explained in Section 2.3. 



The Weyl Ward identities are similarly obtained from Eq. ( 2.20 ), which we rewrite 
here: 

5„W[g] = -2 f ^ 5^^^ (T^,) a = -2 [ dx ^ A{x;g) a{x) , (B.3) 



with anomaly given by ( |A.1| ). Higher-point functions are obtained by differentiating both 
sides of this Equation w.r.t. to g^'^ {J h ■ S/5g): on the left hand side, the variation can 
act on the explicit metric used for the trace or on v.e.v. yielding another stress-tensor; 
on the l.h.s., it can act on ^ or on the anomaly (whose variations are written in the 
notations ( [A.6| , |A.8| )). After two variations, flat-space limit and Fourier transform, one 
finds Eq.(|3.4|), which wc rewrite here: 



{Q{k3)iT{ki).hi){T{k2).h2)) 

+ [{{T{-k2).hi){T{k2).h2)) + {{T{ki).hi){T{-ki).h2))] 

= Aiki,hi;k2,h2)-^[{hi)A{k2,h2) + {h2)A{ki,hi)] . 



(B.4) 



As said in the text, the compatibility of this identity and its analogue for {k^ ^ ki) 
imply the vanishing of the TZ^ term in the anomaly (r = 0), in agreement with the Wess- 



Zumino consistency condition |12|[15|. The compatibility of the Diff and Weyl identities 
yields the equation: 

= A{ki,ki® V + V ® ki]k2,h2) 

+ A{ki+k2,v0(h2.ki) + {h2-ki)^v) 

- {k3.v)A{k2, h2) - {k2.v)A{ki + /c2, /12) - (fei.v)^(/c2, /12) 

- ^ih2)A{ki,v(^ki + ki(g)v) , (B.5) 

which provides an useful check for the variations of the anomaly. The first variation of the 
Euler and Weyl terms vanishes in flat space and thus their second variation is transverse. 
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The conformal isometries of the metric, 

= (,5, + 6^) g^''{x) ^ D^^e^^ - 2a{x)g^''{x) , (B.6) 

imply Ward identities which are hnear combinations of the previous Diff and Weyl iden- 
tities. Let us consider the integrated expressions of the latter equations, say the 0(/i") 
terms: these contain the (n + l)-point function multiplying 6^g^'^ and S^g^'^ , respectively, 
plus the n-point functions and the n-th variation of the anomaly for the Weyl identity. By 
summing the two identities, the (n + l)-point functions cancel out for conformal isome- 
tries ( |B.6| ); after integration by parts of S^h in (B.2) plus 6a-h = —2ah, with vanishing 
surface term |2^, the remaining n-point function contains the "conformal Lie derivative" 
(Je + 5a)T{xi) ( ^.23 ) acting on each stress tensor, z = 1, . . . ,n; finally, there is the term 
with n-th variation of the trace anomaly. The result for n = 2 is given by Eq. ( |3.22| ) . 



C Construction of the basic six-index tensors 

In this Section we construct the tensorial basis on which the correlators are decomposed. 
We start by considering the most general basis of linearly independent tensors with 2n 
indices (n = 1, 2, 3), /^i/i'i • • • fJ-ntJ-'n^ that are symmetric by exchange of any couple Hi <-> 
//'j. In addition we require those tensors to be built in terms of the available covariant 
quantities: A;i^, A;2^, r?^^'. It is evident that all the symmetric 2-tensors made of ki,k2,r] 
are generated by the following iV(2) = 4 2-tensors: 

ihM.h), {k2.h1.k2), {k1.h1.k2), (/ii) (C.l) 

(indices are contracted with the polarization tensor hi as explained in Section 2.3). 

The construction of symmetric 4 and 6 tensors made of ki,k2,r] is simplified by in- 
troducing the idea of primitive tensor: a primitive tensor is a tensor that is made by 
a connected chain of index contractions. For instance {k1.h1.h2.k2) is primitive, while 
{ki.hi.h2.k2){h^) is not. All the 2-tensors of Eq.( p.lD are primitive. 

Here follows the list of the Np{A) = 5 symmetric primitive 4-tensors: 

{k1.h2.h1.k1), {k2.h2.h1.k2), {k1.h2.h1.k2), {k2.h2.h1.k1), (/ii./i2), (C.2) 
and of the Np{6) = 13 symmetric primitive 6-tensors, 

{k1.h3.h1.h2.k1), {k1.h3.h2.h1.k1), {k1.h2.h3.h1.k1), 
{k2.h3.h1.h2.k2), {k2.h3.h2.h1.k2), {k2.h2.h3.h1.k2), 
{k2.h3.h1.h2.k1), {k2.h3.h2.h1.k1), {k2.h2.h3.h1.k1), 
{k1.h3.h1.h2.k2), {k1.h3.h2.h1.k2), {k1.h2.h3.h1.k2), 
{h1.h2.h3). (C.3) 



35 



The complete basis of symmetric 2n-tensors can be constructed by multiplying in all 
possible ways primitive symmetric 2k tensors (1 < k < n) in such a way that the total 
number of free indices is 2n. In particular it is easily checked that the 4-tensors basis 
has dimension iV(4) = 21 = iVp(4) + Np{2f, while the 6-tensors are 137 = Np{6) + 
3 Np{4) Np{2) + Np{2)\ 

For our purposes we also need the basis of 5-tensors -E^iju'^/i2/u'2M3 ™^cle of ki,k2,r], 
that are symmetric by exchange of Hi <-> they are contracted with h\,h2 and one 
polarization vector v. Again, the complete basis is found by multiplying primitive tensors 
in such a way to obtain 5-tensors. In this case, we also need to consider the Np{l) = 2 
primitive 1-tensors: 

{k^.v),{k2.v)- (C.4) 

the iVp(3) = 2 primitive 3-tensors: 

{kiM.v),{k2M.v)\ (C.5) 
and the Np{b) = 4 primitive 5-tensors, 

{v.h2.hi.ki), {k1.h2.h1-v), {v.h2.hi.k2), (fe2./i2-^i-f )• (C.6) 
We easily obtain N{5) = 62 = Np{5)+Np{A)Np{l)+2Np{3)Np{2)+Np{l)Np{2f 5-tensors. 

Next, we need to impose the Bose symmetry realized by the exchange of {ki,hi) <r^ 
(^2,^2) on the 4-tensors and 6-tensors. This can be done by looking at the action on 
the primitive tensors: it is easy to realize then that the Bose symmetric 4-tensors are 
Nf,{4) = 14 and that the Bose symmetric 6-tensors are Nf,{6) = 77. 

In Section 2.3, we defined the tensorial dimension St of a tensor monomial as being 
the total number of ki, k2 with free indices (i.e. contracted with polarization tensors hi or 
v). For instance dT{{kl.h2.k2){k2.hl.k2)) = 4. The tensor dimension of a polynomial is 
the maximal tensorial dimension of the constituent monomials. Hereafter we report the 
dimensions of the tensors basis, divided according to the 6t value, writing m^^, for "m 
tensors of tensorial dimension 6t" ■ 

Classification of primitive tensors: 





= 2 = 


2i , 


Np{2) 


= 4 = 


32 + lo , 


Npi3) 


= 2 = 


2i , 


Npi4) 


= 5 = 


42 + lo , 


Npi5) 


= 4 = 


4i , 


Np{6) 


= 13 = 


= 122 + lo 
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Classification of all tensors: 



N{2) = 4 = 32 + lo , 

Ar(4) = 21 = 94 + IO2 + 2o , 

N{5) =62 = I85 + 323 + 12i , 

N{6) = 137 = 276 + 634 + 422 + 5o . (C.8) 

Classification of all Bose-symmetric tensors: 

iVb(4) = 14 = 64 + 62 + 2o , 

Nb{6) = 77 = 156 + 344 + 242 + 4o . (C.9) 



D Long formulae 



Here we report the explicit form of the "null vector" which occurs in the (1 2) symmetric 
six-index tensors basis Vi, i = 1. ... ,77, (Table |l|) in four dimensions; this vector expresses 
the linear dependence in this basis and should be put to zero (Section 3). The null vector 
is Ts in Eq. (|3.8| ); it reads: 

Ts = 2 (P17 + 7^20-7^26 -7^27 +7^29 -7^33 -7^34 + 7^46 +7^49) -7^36 -7^39 -7^41 -7^43 
+ 2fc2 + V5I + 7^52 - 7^60 - ni - 7^66) 

+ (2/c2 _ g2) + + + + p^_^ + _ _ _ _ _ _ jjgg) 

+ k'^ (-n? - ^8 - 7^64 + 7^69 + 7^70 + 7^72) 

+ ^ (4A;2 - q^) (-2P74 + 7^75 + 7^76 - 7^77) = . (D.l) 

Note that this tensor has 6t = 4, because it possesses no components in the first 15 
elements of the basis; moreover, it is easy to check that its scale dimension is A = 4 



Next, we write the complete set of reduction formulae (Section 4.1) which determine 
the 13 (TTT) amplitudes ^^(g^A;2 = 0) of the St = Q tensors (see Table |): for i = 
1, . . . ,7, they correspond to the traceless transverse tensors, for i = 9, 10, to the anomalies 
(Ae, Ay\r), and for i = 15, 16, 17, to the traceful transverse tensor that vanish at the Weyl- 
invariant critical points. The formulae express the amplitudes in terms of the coefficients 
Ui of the 6t = 6 tensors Vi, i = 1, . . . , 15. 

Three amplitudes and the two constraints among the Uj's were already given in the 
text, see Eqs.( ^ - p^ ); the remaining 7+3 amplitudes are: 



Al = Ui4 , ^2=^2, ^4 = li4 , 
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^3 = TT77 (4^*2 - 27tt3 - 4«4 - U5 + Ouy + llttg - 2nii + 12ui4 + Suis) , 

^5 = ^ (-28u2 + Sl-us + 28ii4 + 7u5 - lAus + buu - lluu - Suis) , 

= -i- {-A0u2 + 99u3 + 40^4 + WU5 - 20us - 18ug + 20un - 12um - 8^15) , 
288 

Aj = m , (D.2) 
and 

Ai5 = - (4n2 - 4m4 - lis + 2m8 + 12^9 + 3nio - 2uii) , 

O 

^16 = ^ (4^2 - 4m4 -U5 + 2us + till + Suis) , 

Aij = — {8u2 — 27n3 — 8^4 — 2^5 + 27^6 + 4^8 — Auu — \2uii — Suis) . (D.3) 

The corresponding formulae for non-vanishing k"^ can be found in the Mathematica note- 
book. At criticahty, there are 3 further constraints ( [4. 81 ) among the Uj's, which are equiv- 
alent to the conditions Ai = 0, for i = 15, 16, 17. 
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